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A general formalism of the problem of perfect state transfer is presented. We show that there 
are infinitely many Hamiltonians which may provide solution to this problem. In a first attempt 
to give a classification of them we investigate their possible forms and the related dynamics during 
the transfer. Finally, we show how the present formalism can be used for the engineering of perfect 
quantum wires of various topologies and coupling configurations. 
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I. INTRODUCTION 

The faithful transfer of a quantum state between two 
distant but specified components of a quantum computer 
is one of the main requirements for practical quantum 
computation T|. The two components (e.g., small quan- 
tum processors) are typically parts of a larger quantum 
network and connected via a quantum channel (wire). 
The state transfer is achieved by converting a stationary 
information carrier (qubit) to a movable ("flying") one at 
the input of the wire [l| . The flying qubit is then trans- 
mitted through the wire towards its other end where it 
is converted back into a stationary qubit. 

In general, photons are excellent flying qubits as they 
can be transmitted coherently over very large distances 
. Hence, many discussions in the field of quantum com- 
puting consider that the state to be transmitted is first 
imprinted onto a photon which is used as the flying qubit 
over an optical fiber. Clearly, the implementation of such 
an idea requires a perfect interface between optical sys- 
tems and the main hardware of the quantum computer. 
However, such an interface is not always an easy task 
to realize as the quantum hardware may be based on 
atoms, ions, molecules or solid state systems, for exam- 
ple. Hence, for short-distance communication (e.g., be- 
tween two quantum processors) it is desirable to develop 
new systems which are fully compatible with the quan- 
tum hardware and, in addition, are suitable for faithful 
state transfer. 

From the point of view of quantum control the problem 
under consideration can be rephrased as follows. Given 
a set of prescribed elements and gates one has to struc- 
ture a system realizing faithful state transfer for arbitrary 
input states. Elements of quantum information pro- 
cessing are typically quantum mechanical objects (e.g., 
ions, quantum dots, Josephson junctions, etc) arranged 
to form linear chains or planar structures with bipartite 
interactions which can be manipulated by certain control 
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parameters. Manipulating the strength of the interaction 
between the quantum elements one can design different 
types of interactions for the whole system and perform 
certain gate sequences. 

Quantum wires based on systems of permanently cou- 
pled quantum objects are of particular interest as they 
require minimal external control thus avoiding signifi- 
cant errors due to the application of multiple operations 
(gates) and/or measurements on various sites of the chan- 
nel. The design of such a kind of passive quantum wires 
and the problem of faithful state transfer have attracted 
considerable interest over the last years. It was imme- 
diately realized, however, that faithful state transfer is 
not an easy task, even in the absence of dissipation or 
dephasing, due to the dispersion of the quantum infor- 
mation along the wire. Thus, various methods have been 
proposed to partially circumvent this problem and to im- 
prove the fidelity of the transfer [3, 4, 5, 6, 7, 8, 9, 10]. 

Throughout this work we focus on perfect state trans- 
fer over passive quantum wires. So far, the existence of 
perfect passive quantum wires has been demonstrated in 
the context of coupled harmonic oscillators llL 12|, ar- 
rays of quantum dots [ISj] and spin chains 10lll4[l5l.[l6t. 
The key idea is the engineering of certain coupling con- 
figurations (and thus interactions) along the chain which 
are able to suppress possible dispersion effects and en- 
force the complete refocus of the transmitted quantum 
information at the ends of the chain at well defined in- 
stants of time. Nevertheless, the majority of these inves- 
tigations mainly focus on centrosymmetric (also known 
as mirror symmetric) linear chains while the Hamilto- 
nian characterizing the quantum wire involves nearest- 
neighbour (NN) interaction only. As a result, they are 
not applicable to realistic setups where these assumptions 
are relaxed (e.g., beyond NN coupfings) Q, and the de- 
sign of new perfect quantum wires is necessary [13]. In 
the worst case scenario where such a design is not fea- 
sible, numerical optimisation techniques may be invoked 
for improvement of the fidelity of the transfer. An alter- 
native solution to this problem relies on the use of the 
so-called dual-rail encoding involving two nearly identical 
quantum channels [l8[. However, in this case the state 
transfer does not occur at well defined time instants and 
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thus the arrival of the state can be revealed by means of 
measurements only. 

The purpose of this paper is manifold. First we aim 
to solve the problem of perfect state transfer (PST) in a 
much broader context. More precisely, we want to inves- 
tigate the possible forms that the system's Hamiltonian 
may take to be suitable for PST. In contrast to previous 
work in the field, here we do not set any a priory restric- 
tions to the topology of the system and the configura- 
tion of couplings between different sites of the channel. 
Hence, irrespective of topology and coupling configura- 
tion, we show that there are infinitely many Hamiltonians 
which are suitable for PST. For the sake of illustration, 
we demonstrate that certain Hamiltonians previously dis- 
cussed in the literature can be obtained in the framework 
of our unified theory by setting certain restrictions and 
using the right parameterisation. Second, we show that 
our approach provides new ways for quantum wire en- 
gineering for systems of arbitrary topology and interac- 
tions beyond nearest neighbours, thus generalising exist- 
ing work in this context [3,[H)[l3- However, in contrast 
to [H, [H, our approach to the problem of PST docs 
not rely on the concept of inverse eigenvalue problems, 
but rather on the derivation of Hamiltonians which lead 
to a particular transform namely, a permutation opera- 
tion. Finally, the present work can also be viewed as an 
attempt to classify the Hamiltonian problems leading to 
exact revivals in discrete time evolution and the types of 
dynamics that appear during this evolution. 

Our theoretical approach was also used for the deriva- 
tion of Hamiltonians for the Fourier transform imple- 
mented using linear optical elements [l^. The general 
form of the Hamiltonian was given and the results are 
closely linked to the present ones as the Fourier trans- 
form is a cyclic operation (a property shown to be crucial 
for our considerations) . The square of the Fourier trans- 
form gives a permutation operation, hence is related to 
the problem under investigation in the following sections. 

The paper is organised as follows. In Sec. |TT]we intro- 
duce our mathematical model and derive a class of Hamil- 
tonians suitable for PST. In Sec. IIIIl considering small 
networks we demonstrate how our theoretical approach 
can be used for quantum wire engineering. We show that 
certain Hamiltonians previously discussed in the litera- 
ture are members of the larger unifying class derived in 
Sec. ini Moreover, we derive new PST Hamiltonians in 
the framework of NN-type, coulomb and dipole-dipole 
interactions. In Sec. llVi we analyse and discuss the dy- 
namics induced by various types of PST Hamiltonians 
while we conclude with a summary of our main results 
in Sec. B 



II. HAMILTONIANS FOR PERFECT STATE 
TRANSFER 

The model we consider is a network consisting of n 
sites labelled by {1, 2, . . . , n} = S„. At time t = 0, the 



qubit (e.g, spin) in the 1-st (input) site of the network is 
prepared in the state IV'in)- We wish to transfer to the n- 
th (output) site of the network with unit efhciency after 
a well defined period of time, let us say t-c = r/ J, where 
J^^ is our time units (fi. = 1). 

We are interested in perfect passive quantum wires 
i.e., networks involving permanently coupled sites with- 
out any additional external control. Depending on the 
physical realization of the network, J can be a char- 
acteristic energy in our system, a coupling constant or 
tunnelling (hopping) rate, etc. Nevertheless, to present 
a generic theoretical approach to the problem of PST, 
the following discussion will be based on the dimension- 
less quantity r. It is also worth noting, that r need not 
always be a continuous variable as in some physical real- 
izations (e.g., passive linear optical networks), the excita- 
tion (i.e., photon) evolves under successive applications of 
identical (practically instantaneous) unitary operations 
induced by the system's Hamiltonian [20]. In this case 
T is simply the number of unitary operations we have 
applied. The subsequent discussion applies to all these 
physical realizations no matter whether they involve dis- 
crete or continuous evolution. However, for the sake of 
brevity we mainly refer to r as number of applications 
(iterations) of identical unitary operations. 

Our purpose is to explore the possible forms that the 
system's Hamiltonian Ti. may take, to be suitable for 
PST. In other words, we wish to analyse the whole class 
of Hamiltonians which lead to PST from the first to 
the last site after exa ctly r applications. Following previ- 
ous work in the field [ifl [H H E Q [H E [III, III , we 
will assume that the network does not disturb the trans- 
mitted state and the main source of possible dissipation 
and decoherence is the fact that the wave packet (quan- 
tum information) spreads along the network. Hence, we 
are interested in Hamiltonians which preserve the total 
number of excitations in the system. For instance, in 
the framework of spin chains and arrays of quantum dots 
this means that the Hamiltonian commutes with the total 
spin operator [H [H, [H [H, [13, 111 . 

One way to guarantee all of these requirements is to 
concentrate our subsequent investigation on Hamiltoni- 
ans for which the associated unitary evolution lA, leads 
to a permutation matrix after t applications i.e., 



U{t) 



where 



V = 



/O 



\i ... oy 



(1) 



(2) 



is a permutation in the single excitation sector i.e. also 
the submatrix V in the computational basis is a per- 
mutation. The permutation matrix ^ guarantees that 
the net effect of the evolution ([T]) is the transfer of the 
excitation from the first position to the last one i.e.. 
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|l)i|0)2...|0)„ ^ |0)i|0)2...|l)„. The condition ^ might 
seem restrictive at this point and this is in general indeed 
the case as one may derive Hamiltonians which satisiy the 
aforementioned requirements and do not lead to permu- 
tations, but rather to other unitary operations. Never- 
theless, as we show later on, many known Hamiltonians 
suitable for PST are basically associated with permuta- 
tions and can thus be obtained within the present unify- 
ing theoretical framework. Besides, we will demonstrate 
that we have infinitely many, yet unexplored, choices for 
quantum wire engineering which are not covered by pre- 
vious work in the field. We have to point out, however, 
that condition ((!]) is indeed a severe restriction if one is 
interested only in perfect transfer of a particular set of 
possible states or properties of states such as probabili- 
ties, entanglement, etc. 

Note now that perfect transfer of arbitrary single-qubit 
states is a sufficient condition for perfect transfer of ar- 
bitrary multi-qubit states. For instance, in this case the 
same perfect quantum wire can be used multiple times 
for the transfer of each qubit separately. Furthermore, 
as long as the number of excitations and the transmitted 
qubit state are preserved by the Hamiltonian, the total 
Hilbert space in the problem can be decomposed into sub- 
spaces and we can focus on the one-excitation subspace. 
The corresponding basis is denoted by {\a) | a S 
and indicates the presence of the excitation at the site 
a. Hence, the problem of the PST is essentially reduced 
to the perfect transfer of a single excitation initially lo- 
cated at the 1-st site i.e., |V'in) = 

In general, we can define (n— 1)! different permutations 
of the form we have n — 1 free positions to fill, and 
for each one of these our purpose now is to construct 
the class of Hamiltonians satisfying Eq. ([T]) and thus are 
suitable for PST. We can distinguish between cases where 
V consists of one or more cycles. In the following, we are 
going to investigate the two cases separately but we will 
see that the latter reduces to the former one within each 
cycle. 



A. One-cycle permutations 

For n sites we have (n — 2)! possible one-cycle permuta- 
tions. In particular, starting from the n x n permutation 

/O 1 ••• 0\ 



\ 1 0/ 

the other possible one-cycle permutations can be ob- 
tained by simply relabelling the (r7,~2) intermediate sites 
{2,...,n— 1}. In the following we will derive our results 
for this particular matrix. 

To determine the class of Hamiltonians C-^ satisfying 
Eq. ((T|) for the permutation ([3]), we need to know the 



eigenvalues and eigenvectors of the operator V. The spec- 
trum (T consists of n different eigenvalues of the charac- 
teristic equation A" = 1 i.e., 

a ^ {Xj \ Xj ioi- i,j eZn}, (4) 

where 

A," = exp ( i27r— ) for j e Z„, (5) 
V "/ 

and Z„ = {0, 1, . . . , n— 1}. The corresponding normalised 
eigenvectors can be expanded in the computational basis 
{\a) I a G §„} as follows 

and thus the spectral representation of V is 

A J go- 
Using Eqs. (I5|)-(l7|), it is straightforward to construct a 
first Hamiltonian which satisfies Eq. ([T]) and ([3]) as follows 

^ arg(A,)|yA,)(2;A,| (8) 

Aj Go- 
where arg(Aj-) is the phase of the j-th eigenvalue Xj. 
However, this is not the only Hamiltonian which leads 
to the permutation V after r applications. Additionally, 
wG may shift each eigenenergy of the Hamiltonian ([8]) by 
an arbitrary integer multiple of 27r obtaining 

Aj go- 

= 7 E MA,) + 2^?^^] |2/A,)<yA,|, (9) 
with 1 e Z" = {(/ao,;ai,...,/a„_i) I e Z,j e Z„}. 

Therefore the class of Hamiltonians satisfying Eq. ([T|) for 
the particular one-cycle permutation ([3]) is 

= {7^1 I 1 e Z"}, (10) 

with Til given by Eq. ([5]). Since the members of <Lh are 
parameterised by the integer vector 1 e Z", the class 
<Lu consists of infinitely many Hamiltonians suitable for 
PST. 

Note now that all the members of (Z-u are linear su- 
perpositions of the projectors onto the eigensubspaces of 
the permutation V. More precisely, Eq. ([9]) can be also 
written as 

He = X! \y^^) (y^^= E £^.nA,' (11) 

Ajgo XjG(T 

where the spectrum (eigenenergy-vector) E G R" = 
{(eao, ■ ■ • ,£a„_i) I £\, e R,j £ Z„}. 
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Using Eq. ([6]), we can express Eq. (|lip in the compu- 
tational basis as follows 

E T.^>^.^r'\^)^p\^ (12) 

a,/3eS„ Xj So- 

which can be also written in the usual form 

We = We ^ + Ve (13) 
with the diagonal part 

n'i^ = E I") ("I (14) 

and the interaction 

Ve- E G(a,/3)|a)(/3|. (15) 

Thereby, the energies and the couplings are given by 

= E (16) 

G(a,/3) = E^^.^""''- (17) 

From Eq. ([TB)) . we see that Hamiltonians which are asso- 
ciated with permutation ([3]) and lead to PST may only 
correspond to networks involving the same energy level 
for all the sites. 

The conditions under which Hamiltonians with 
NN interaction may lead to PST have been exten- 
sively discussed in the literature during the last years 
[S El [H, Q m, [Hi- Such NN-interaction Hamil- 
tonians are typically tridiagonal in the computational 
basis i.e., we have G{a,(5) = 0, for all (5 ^ {a, a ± 1}. 
However, according to the following theorem the class of 
Hamiltonians £-h we have derived here does not include 
NN-interaction Hamiltonians. 



where Vx involves NN interaction i.e., /3 = a± 1, while 
Vx'^^ involves terms with (3 ^ {a, a± 1}. Our assumption 
about the existence of an NN-type Hamiltonian implies 
that the Hamiltonian is tridiagonal in the computational 
basis, i.e., 

V^^) ^ and V4^' = 0. (19) 

From the second condition and Eq. (|15p we have that 
G{a,f3) — 0, for all /3 ^ {a, a ± 1}. However, since 
G(a, P) depend only on a — /3 we may obtain n — 2 linear 
equations for variables Xj (see Eq. [T7|) . namely Ax = 
where 

/ Aq Xl ■ ■ ■ Xn^i \ 
A= . . . (20) 

\n— 1 \n— 1 \n— 1 

\ -^0 -^1 ■ ■ ■ I 

Note now that, in view of Eq. ([5]), the matrix (PO)) is 
(up to first two missing rows and normalisation) the 
usual Fourier transformation. Hence, rank(A) = n — 2, 
while two linear independent solutions can be chosen as 
xi — (1, . . . , 1) and X2 (Ao, • • ■ , A„_i). Asking for a 
linear combination of xi and X2 to be real (because this 
plays the role of eigenenergy vector) we find that the 
only acceptable solution is multiple of (1, . . . , 1). How- 
ever, one can see immediately from Eqs. and 
that such a spectrum can be associated only with the 
system involving no interaction between sites i.e., 

Vf^ = and ' = 0. (21) 

which contradicts our initial assumption (fT9|) about NN 
interaction. □ 

Closing the section we would like to emphasise that the 
presented method (including the proof) applies also to all 
other one cycle permutations. The explicit results will 
differ from the presented ones only by the corresponding 
permutation of the labels. 



Theorem. For networks of arbitrary dimension 
{n > 2), there exists no nearest-neighbour-interaction 
Hamiltonian satisfying condition {Ip in the framework 
of permutation Q)- 

Proof. According to Eq. (fTT|) all the Hamiltoni- 
ans which satisfy Eq. ([T]) with V given by Eq. ([3]) are 
linear superpositions of the projectors Hj. Moreover, 
they can be decomposed as described in Eqs. P^ - P?)) . 
For a given dimension n > 2, let us now assume that 
there exists a set of eigenenergies x £ R" such that the 
corresponding Hamiltonian Tix involves NN interaction. 
We are going to show that such a Hamiltonian does not 
exist (proof by contradiction) . 

We start with the general decomposition of the inter- 
action term Vx as 



Vx = v(^) + v(^) 



(18) 



B. Many-cycle permutations 

In the case of permutations consisting of more than 
one cycle, the situation is slightly more complicated but 
the problem can be treated separately within each cycle 
along the lines of the previous section. More precisely, 
consider a cycle of length d < n. Such a cycle does not 
involve all the network sites but rather a subset of them 
Sd C §„. Hence, the related eigenvalues are given by 



Ao ~ exp I i27r— 
* d 



for j eZd, 



(22) 



while for the projections of the corresponding eigenvec- 
tors onto the computational basis {\a) | a G §„} we have 



Xr^/Vd for aeSd 



otherwise 



(23) 
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As is apparent from Eq. ([22)1 the spectrum tr of a many- 
cycle permutation is always degenerate since at least one 
of the eigenvalues (i.e., the eigenvalue for j = 0) appears 
as many times as the total number of cycles. Hence, the 
same degenerate eigenvalue Xj corresponds to Sx . dis- 
tinct eigenvectors from different cycles. To this end, in 
Eq. ((23|) the eigenvectors are characterised by the ad- 
ditional superscript k € {1, . . . , }, where is the 
degeneracy of the eigenvalue Xj. For a given degener- 
ate eigenvalue Xj , let us also denote by the subspace 

spanned by the S\- distinct eigenvectors {l^'i^^)}- 

In view of the degeneracy in the spectrum of a many- 
cycle permutation we have more freedom in the construc- 
tion of Hamiltonians satisfying condition ([T]). Indeed, as 
described in Sec. Ill A[ for such a construction we first of 
all need a basis of orthonormal eigenvectors. However, 
within each eigensubspace £\. we can choose such an 

eigenbasis {ly^'^!'')} in many different ways, by construct- 
ing linear superpositions of Moreover, for each 

eigenvector \y'"x^) we have the additional freedom to shift 
the phase of its eigenvalue by a multiple of 27r, as we did 
in Sec. Ill Al So, the class of PST Hamiltonians ([10]) has 
infinitely many members of the form 



Xj k=l 



(k) 



y^) 



Aj k=l 



(k)\ / (fc) 

y\;){y\; 



(24) 



where leZ"^{(^«,...,/i^^ 

Z}. This is a straightforward generalisation of Eq. ^ 
to the case of degenerate eigenvalues. 



C. Quantum wire engineering 

So far we have seen that given a permutation matrix 
of a particular form one may easily derive the whole class 
of PST Hamiltonians which lead to the permutation un- 
der consideration. Hence, our method provides a way 
of quantum wire engineering in the sense that it enables 
us to construct possible PST Hamiltonians (i.e., to de- 
fine energies and coupling strengths) which are imple- 
mentable in a particular setup. The starting point is 
always the choice of the permutation which need not be 
always random, as various factors (such as topology of 
our network and type of interactions) may automatically 
fix our permutation matrix. For instance, in the case of 
a centrosymmetric network one has to focus on antidiag- 
onal permutation matrices only while, according to the 
theorem of Sec. Ill Al NN-type Hamiltonians automati- 
cally exclude all one-cycle permutations. 

Having fixed our permutation, the construction of a 
PST Hamiltonian proceeds in two steps. In the first step 



one has to define the corresponding class of PST Hamil- 
tonians working along the lines of this section. The mem- 
bers of this class are parameterised by a number of free 
parameters. In the second step one can estimate all these 
parameters by applying certain constraints based on the 
topology of the network as well as the form of the physi- 
cal interactions implementable within the framework of a 
particular setup. This engineering process can be always 
performed numerically but in certain cases (especially for 
relatively small networks) derivation of analytic solutions 
might be possible. In the following section we discuss 
the quantum wire engineering in detail by explicitly con- 
structing PST Hamiltonians for small networks. 



III. EXAMPLES OF PST HAMILTONIANS 

Let us devote some attention to the application of the 
proposed quantum wire engineering method to few con- 
crete examples. Among the simplest we can think of, is 
the example of a small networks consisting of just few 
sites. 



A. Nearest-neighbour interaction 

Following the notation introduced in the previous sec- 
tion, we start with the investigation of the whole class of 
NN-type PST Hamiltonians in the context of a small net- 
work consisting of four sites only. The problem of PST 
in such a small network is amenable to analytic solutions 
thus offering the appropriate theoretical framework to 
demonstrate the application of our method. 

We can define six 4x4 permutation matrices of the 
form (m that is, two one-cycle permutations and four 
permutations involving more than one cycles. As we 
showed in the previous section, NN-type Hamiltonians 
cannot be obtained in the context of one-cycle permuta- 
tions. Hence, for the purposes of this particular example, 
it is sufficient to restrict ourselves to permutation matri- 
ces involving more than one cycles only. 

Let us consider first the antidiagonal permutation ma- 
trix 



1 
10 
10 
10 



(25) 



This permutation involves two cycles of dimension d = 2 
namely, Ci = (1,4) and C2 = (2,3) [l^. Hence, accord- 



ing to Eq. ([22)1 we have two eigenvalues Aq = e = -1-1 
and Al = e"^ = —1. Both eigenvalues ±1 are dou- 
bly degenerate (i.e., S± — 2) and the corresponding 

eigenvectors are given by |w±^) = (|1) i |4))/V2 and 

[■^i^) = ± |3))/-\/2, for cycles Ci andC2, respectively. 
Accordingly, the corresponding subspaces are 



,(1) 



,(2) 
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Having defined the permutation matrix, the construction 
of PST Hamiltonians now proceeds in two steps. 



1. Step 1 — Parameterisation. 

For each one of the subspaces we can choose an or- 
thonormal basis in many different ways. For example, 
taking into account the orthogonahty condition, we may 
choose for the two subspaces 



(1) 



„(i) 



,(2) 



and 



,(2) 



(2) 

y- 



,(1) 



,(2) 



(26) 



,(2) 



,(1) 



,(2) 



Thereby, ^, i^, ^, C G C such that 



1 and |^|2 + |C|2^f. 



(27) 



(28) 



According to the theory of the previous section, the cor- 
responding eigenergies can be chosen as 



+ 27r/^^^ 
+ 2ttI 



+ 

(2) 



(29) 



and 



-(1) 



TT + 2t:1 



(1) 



(2) 

e_ = TT 



2Trl 



(2) 



(30) 



(k) 

for e Z. Thereby, note that the eigenvalues of sub- 
space 5+ a-re even (odd) integer multiples of tt. 
Thus, the whole class of PST Hamiltonians ([M)) reads 



1 



-(1) 



-(2) 



(2)\ / (2) 

y+ ){y+ 



-(2) 



(2) 

y- 



(2) 

y- 



,(31) 



with open parameters the spectrum {e^^* , e*^'', e'^\ e'^'' } 

[or equivalently the integers 1^.' ] and two independent 
complex numbers (due to normalisation), say and ^. As 
we discuss in the following subsection, all these parame- 
ters can be specified by imposing additional constraints 
on the form of the resulting Hamiltonian. Moreover, we 
would like to note that, due to the form of the eigenvec- 
tors, all the members of this class are symmetric along 
the diagonal and the antidiagonal. 



2. Step 2 — Parameter estimation. 

Assume now that we are interested in the whole class 
of NN-type Hamiltonians leading to the antidiagonal per- 
mutation matrix (j25p and are of the form 



7^(4) 



/Si 51 

51 E2 52 

52 E2 gi 

V 51 Si 



(32) 



Since the class ([51]) has been derived within a rather gen- 
eral framework (limited only by the particular form of 
the permutation), our task reduces to the application of 
specific constraints on the parameters entering 7i|"''' . 

We will focus on the case of Hamiltonians with 

(1 2) 

non-degenerate spectrum i.e., all the eigenenergies e^.' 
are different. The case of degenerate spectrum can be 
treated similarly and leads to NN-type Hamiltonians 
with vanishing couplings. Hence, such solutions corre- 
spond to a broken network and can never lead to PST 
from the first to the last site. 

Rewriting Eq. pip in the computational basis and ask- 
ing for Hamiltonians of the form ((32l) . we have that the 
matrix elements (l|7i|^'' |3) and (l|7i|'^-' |4) must vanish. 
Thus we obtain the following set of nontrivial constraints 
on the free parameters, 



e^+^U4'-s^-'\e-e'^'\C\-'^0, (33a) 



-(1) 



(2), 



(33b) 



Other equations which can be obtained from the remain- 
ing matrix elements {i \ 7i|^^ |j) will be used later on, for 
the derivation of sets {Ei, E2} and {51, 52} for which the 
Hamiltonian ^ enables PST. 

Using the normalisation conditions ([28]) , Eq. (|33aP can 
be also expressed in terms of differences of eigenvalues 
only i.e.. 



(1) 



-(2) 



(1) 



.(2) 



-(2) (2) 



0. 



This fact indicates that the existence of a solution to 
the problem of PST depends mainly on differences of 
eigenenergies. In other words, if we have a solution for 
a particular choice of the spectrum, say {si, £2, £3, £4}, 
then there also exist the solutions for all the spectra of 
the form {ei-|-27rj, e2 + 27rj, £3-f27rj, £4-|-27rj} with j e Z. 

Assume now that i' and ^ are real and non-negative. 
Such an assumption is not restrictive because the global 
phase of each eigenvector can be chosen arbitrarily with- 
out changing the Hamiltonian (and thus dynamics). The 
phases of /i and C will be denoted by ip and respec- 
tively. Solving Eqs. p3p with the normalisation condi- 
tions (|28p. we obtain for the absolute values of the pa- 
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rameters ^, and C 



1^1 



ICI 

















)(eV -£ 
















i'l ^ 






- eL^ 4 


-ef) 


















- eL^) 4 




















- 4 





(34) 



Moreover, from Eq. p3bp we have that a solution exists 
if and only if the phases of the amplitudes satisfy the 
relation 



if — X — TTiTT, m e Z 



(35) 



where m must be odd(even) for (e^^'^ — e+ ■*)/(£ 
negative(positive), respectively. Finally, it is straightfor- 
ward to show that the normalisation condition and the 
non-negativity of the right-hand sides in Eqs. 



-(1) 



■ £_ ) 



-(1) J2) 



- + 



and 



-(1) ^(2) 



ply that the intervals between 

must overlap but need not be one inside the other. For in- 
stance, one can readily check that the spectra discussed in 
many papers [H, 0, [l^, [l^ fulfil this condition and thus 
can be obtained within the present formalism. Wc have 
also verified this fact numerically for larger networks. 

At this point, we have analysed completely the PST 
Hamiltonians which lead to the permutation matrix ()25|) . 
Indeed, using the expressions for the absolute values and 
choosing the phases and the spectra according to the 
aforementioned constraints, it is straightforward to de- 
rive possible combinations of {Ei, E2} and {gi, 32} which 
lead to PST transfer by means of the equations 



(4) 



and 



1|hP^|2)=5i, 



mni'^ |3) 



.92, 



respectively. 

The remaining question is whether there are other 4x4 
many-cycle permutation matrices generated by NN-type 
Hamiltonians. To answer, we have to analyse the three 
remaining configurations of cycles. 

• Permutation (4, 1)(2)(3). In this case we have two 
cycles of lengths two, one, and one, respectively and 
two different eigenvalues namely, 4-1 (degeneracy 
S+i = 3) and —1. Following the theory of Sec. IIIBI 
and working as before for the subspace £+, we can 



define new eigenvalues e^,!'' = 2t:1'^' with the cor- 
responding eigenvectors in the computational basis 



7W 



given by \y^^^) = + \4))/V2 + |2) 4- 6 |3), 

for i e {1,2,3}. In the subspace f_, the eigen- 
value and the eigenvector remain unchanged i.e., 
£_ = -1 and |?/_) = (|1) - |4))/^/2. Hence, the 
general form of the Hamiltonian which leads to the 
permutation under consideration after r applica- 
tions is given by 



n 



(4) 



■ 3 

.1=1 



+ 



\v-) {y- 



Recall now that we are interested only in NN-type 
Hamiltonians i.e., Hamiltonians with (4|7Y|'*-' |2) ~ 
and (1|7Y|^^ |3) = 0. Rewriting Ti^f^ in the com- 
putational basis, we find that these two conditions 
also imply {l\'H^f' |2) = and (4|7Y|''' |3> ^ 0. In 
other words, the first and the last sites are com- 
pletely decoupled from their respective neighbours 
and thus there are no NN-type Hamiltonians which 
lead to PST in the context of the permutation 
(4, 1)(2)(3). Strictly speaking, PST from the first 
to the last site may occur only if these two sites 
are directly coupled. This is, however, the case of 
a trivial two-site network for which PST is always 
possible. 

• Permutation (4, 1,3) (2). In this case we have two 
cycles of lengths three and one, respectively. The 
general form of eigenvectors in the computational 
basis is \y^^^) = + |3) + |4))/V3 + |2) and 

\yf) = -^^*(|l) + |3) + |4))/^/3+^i* |2) for the eigen- 
value -hi, \yx) = (|1) + a: |3) + |4))/V3 for the 
eigenvalue x, and = (|1) 4- x"^ |3) 4- x |4))/\/3 
for the eigenvalue a:^, where x — exp(2i7r/3). Ask- 
ing for a NN-type Hamiltonian, among others, we 
obtain the constraints 



-(2), 
-(2): 



£ I £ ^2 X — 



which should be satisfied simultaneously. This im- 
plies the condition {ex — £x^){x ~ x'^) — which 
cannot be fulfiled because Sx and Ex^ are of the 
form Sx = (2/3 4- 213)tt, and £^2 = (4/3 -I- 214)tt, 
with l^jli e Z. 

• Permutation (4, 1, 2) (3). This case is similar to the 
permutation (4, 1,3)(2). 

Hence, for the four-site network, NN-type PST Hamilto- 
nians involving permutation transformation can be ob- 
tained only in the framework of the antidiagonal permu- 
tation matrix ((25|) . 

Concluding this simple example it is worth keeping in 
mind the aforementioned observation about the overlap 
of the intervals for odd and even eigenvalues. This turns 
out to be a particularly useful result as it enables us to 
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construct more PST Hamiltonians by changing appropri- 
ately the spectrum of an ah'eady known PST Hamilto- 
nian. Most importantly, we have verified numerically the 
functionality of this idea for interactions beyond nearest 
neighbours and for networks involving more than four 
sites (e.g., see TableUin the next example). We turn now 
to discuss PST Hamiltonians beyond NN interaction. 

B. Beyond nearest-neighbour interaction 

In the previous example we explicitly discussed the 
construction of NN-type PST Hamiltonians. In the 
problem of PST Hamiltonians beyond NN interaction 



was addressed for the first time. Motivated by this work 
we have studied the problem of PST Hamiltonians with 
a prescribed drop off of the coupling constants also in 
our theoretical framework. In this case, naturally, we 
do have only a very limited control over the form of the 
Hamiltonian. However, before we proceed to present re- 
lated results, we point out that the class of Hamiltonians 
derived in Sec. Ill Al is already of the beyond nearest 
neighbour type. 

We have applied our method to Hamiltonians involving 
two different types of interactions and for networks with 
a moderate number of sites namely, four and six. For the 
latter case we present some of the obtained numerical 
results. Consider the Hamiltonian of the form 



7^(6) ^ 



El 
1 

7 



(ri+ra)^ 
1 

(ri+r2+r3)'< 

1 

{ri+2r'i+r3)-< 



E2 
J_ 

r"' 

i 

{r2+r3)~< 
1 



1 1 1 1 \ 

(ri+r2Y< (ri+r2+r3)i (ri+2r2+r3)^ (2ri-|-2r2+r3)^ \ 
1 1 1 1 



1 

7 



V . . - ^ . • ... ■,- . . - ^j (r2+r3)-' 
\ 1 1 1 

\ (2ri+2r2+r3)T (r-i+2r2-|-r3)T- (ri+r2+r3)"' 



(1-2 +1-3)"' 
1 



Ez 
1 

(ri+r2)"' 



(2r2+r3)-i 
1 

(r2-|-r3)^ 
1 



E2 
1 



(ri-|-2r2+r3)T 
1 

1 

{ri+r2)~< 
1 



El 



which refers to a network where interactions extend be- 
yond nearest neighbours and drop off with distance as 
l/r'*'. Our purpose is to find combinations of energies 
i?2, ^^3} and distances {ri,r2,r3}, such that "H^^^ 
leads to an antidiagonal 6x6 permutation matrix after 
precisely r applications. 

As in the previous example, the Hamiltonian design 
proceeds in two steps. First we have to find a general pa- 
rameterisation of the possible PST Hamiltonians along 
the lines of Sec. IIIBI Subsequently, asking for Hamil- 
tonians of the form p6p . we obtain a set of equations 
for the free parameters entering our design. This sys- 
tem of equations can be solved by means of standard 
numerical techniques yielding a number of acceptable so- 
lutions (i.e., sets of energies {£^1, £^2, £'3} and distances 
{ri, r2, r^}). A few of them, for the coulomb (7 = 1) and 
the dipole-dipole (7 = 3) interaction are given in Ta- 
ble HI In discussing this table, we would like to note that 
the solution (a) for the dipole-dipole interaction coincides 
(apart from a shift in the energies) with the one given by 
Kay while all the other solutions are new and have 
not been discussed in the literature before. Moreover, 
as is evident from the table footnotes, for both solutions 
(a) and (6) we have an overlap of the spectrum for even 
and odd eigenvalues. More solutions can be obtained 
in a straightforward manner by modifying appropriately 
these spectra, while they overlap. 

Throughout this section, for the sake of simplicity and 
illustration, we have focused on small networks and par- 
ticular choices of Hamiltonians. In closing we would 



TABLE I: Perfect state transfer. Parameters entering the 
Hamiltonian (|36|l , for 7 = 1 (Coulomb interaction) and 7 = 3 
(dipole-dipole interaction). The depicted values for the en- 
ergies and the distances should be multiplied by vr/r and 
(t/tiY^~' , respectively. 





Coulomb interaction 


Dipole interaction 


Parameter 


Solution" 


Solution' 


Solution" 


Solution' 


El 


1.07571 


1.10848 


-0.00885 


0.36328 


E2 


-0.70069 


-0.73164 


-0.61799 


-1.89782 


E3 


-1.87501 


-3.87684 


-0.87315 


-1.96546 


ri 


1.77185 


1.93029 


0.96704 


0.97927 


r2 


1.12358 


0.71338 


0.90156 


0.74773 


r-i 


0.96133 


1.06286 


0.88587 


0.89090 



"Eigenenergies: {si , £2, £3, £4, £5, £5} = {2, 0, —2, —1, —3, 1}. 
'Eigenenergies: {si, £2, £3, £4, £5i£6} = {2, 0, —4, —1, -5, 1}. 



like to emphasise that the presented solutions are by 
no means exhaustive. Our theoretical approach to the 
problem of PST is able to provide infinitely many solu- 
tions and is not a priori restricted to particular topologies 
and/or coupling configurations. In the following section 
we turn to discuss the evolution of the excitation under 
the influence of various types of PST Hamiltonians. 
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IV. WAVE PACKET DYNAMICS 

The Hamihonians given in the previous sections are 
quite abstract. To illustrate in more detail the dynamics 
they can induce we present a few examples of how the 
dynamics of a single excitation will look like. In par- 
ticular we will comment on two parameters namely the 
occupation probability distribution and the correspond- 
ing entanglement. As we mentioned before, there are 
infinitely many PST Hamiltonians one may propose, and 
each one of them may induce a new type of dynamics. 
Hence, the patterns we have chosen to discuss in this 
section at any rate cannot be considered as representa- 
tives of the whole range of patterns generated by all the 
possible PST Hamiltonians. 



A. Occupation probability distribution 

In Sec. |TT] we have explicitly constructed the whole 
class of Hamiltonians €t-c which satisfy condition ^ . In 
particular, we have seen that there are infinitely many 
Hamiltonians Tii which in general may lead to different 
evolution operators and thus to different dynamics. A 
more detailed classification is not an easy task. Never- 
theless, in this section we discuss typical transfer dynam- 
ics associated with particular choices of the spectrum, or 
equivalently the integer vector 1 entering the members of 

An excitation initially localized at site i evolves under 
m applications of the Hamiltonian according to 

Ui{m)\i) = e'^""|i) 



Aj Go- fc=0 



(37) 



The corresponding probability to find the excitation at 
site / is given by 



Aj Gcr fc=0 



(38) 



Hence, during the transfer the probability of finding the 
excitation may spread among all the sites of the network 
which become essentially entangled. However, no matter 
how complicated the transfer dynamics may be, we have 
the desired perfect transfer of the excitation after exactly 
T applications since U\(t) = V . 

Let us have a more detailed look at the single exci- 
tation (e.g., photon) wave packet dynamics, initially lo- 
calised at the first site i.e., \i) = |1) in Eq. In Fig. 
[T]we plot the evolution of the occupation probability dis- 
tribution as a function of the number of applications m. 



for a medium size network of 11 sites (modes). The four 
plots differ by the choice of the spectrum. More pre- 
cisely, the first plot has been obtained for the spectrum 
with Zaj = ll^ii the second plot is for the spectrum with 
^Aaj-i = 19 — 3j and = (for 1 < j < 6), and the 
third plot corresponds to the symmetric spectrum with 
= li — 6|. The last plot corresponds to an alternating 
sequence of eigenenergies with /a4j_3 = = 0, and 

^A4,_i = ^A4j = 5, for 1 < j < 3. 

The descendent sequence of eigenvalues used in Fig. 
[ija) results in a gradual transfer of the excitation from 
the initial position to the target position. At each itera- 
tion instant the wave packet is slightly shifted and only 
few of the adjacent modes is noticeable excited. On the 
contrary, for all the other plots the wave packet is not 
only broadened and propagating, but splits into two or 
more components. As a result there is a nonzero proba- 
bility to detect the excitation at two or more distant sites 
at specific instances of the evolution. In particular, when 
the decreasing sequence of eigenvalues is interrupted by 
zero values, the probability distribution pattern is formed 
by several "paths" which transfer the initial excitation to 
the final target site [see Fig. [Ijb)]. When the eigenvalue 
sequence exhibits a symmetric dip, the probability dis- 
tribution is formed by two intersecting lines along which 
the probability propagates [see Fig. [Ijc)]. The final plot 
[Fig. md)], applicable for the tooth like sequence of eigen- 
values, shows an effective transfer of probability along 
the outer parts of the network. The central modes [23| 
are never significantly populated during the wave packet 
propagation. 



B. Entanglement 

The behaviour of the single excitation (e.g., photon) in 
a network implementing PST can be discussed also from 
the point of view of entanglement. We can ask how the 
propagation of the photon in different types of networks 
affects the mutual (bipartite) entanglement. The bipar- 
tite entanglement can be quantified using the concept of 
concurrence [2l|. For the single photon excitation the 
bipartite entanglement between two selected modes, say 
i and j, is given by [l^l 



This simple result shows that shaping the probability 
distribution strongly infiuences the entanglement prop- 
erties. The more the photon spreads across the modes 
during the evolution the more entangled the modes be- 
come. When analysing the bipartite entanglement be- 
tween the initial mode and the remaining modes we find 
that the entanglement plots exhibit a very similar struc- 
ture to the one of the probability distribution shown in 
Figs. [Il^a-d). The main difference is that, due to the 
probability multiplication in the definition of the concur- 
rence, the peaks become broader and less pronounced. 
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Due to the similarity we do not present any plots regard- 
ing the bipartite entanglement. Instead, we proceed to 
discuss the single-excitation dynamics from the point of 
view of the overall bipartite entanglement T(m), where 
additional conclusions can be drawn. 

By definition, T(m) is the sum of all Cfj and is given 



T{m) = 2 



Therefore, we see that the overall bipartite entanglement 
after m applications of the Hamiltonian is specified by 
the squared probabilities (essentially the purity - lin- 
earised entropy related to the probability distribution). 
The more localized the distribution is, the less overall bi- 
partite entanglement the propagating wave pattern con- 
tains. In particular for the patterns presented in Figs, 
[ija-d) we come to the following conclusions concerning 
the overall bipartite entanglement. 

The choice of the spectra influences the maximum 
value of T which can be attained during the propaga- 
tion as well as its modulation. As is depicted in Fig. [2l 
the overall bipartite entanglement for the monotonously 
decreasing [Fig. E^a)] and the symmetric-dip spectrum 
[Fig. HJc)] exhibits simple forms. In both cases the exci- 
tation occupies at least one of the intermediate sites with 
probability one [see Figs. [Ha), (c) respectively] and this 
fact is reflected in T which drops to zero. 

The overall bipartite entanglement for the two other 
spectra acquires a more complicated structure and the 
maximum value of T in Fig. [2Ib) comes rather close to 
the maximum achievable value of 1.81 This just 

indicates that the wave packet is quite uniformly dis- 
tributed among the sites [see Fig. [Hb),(d)]. The question 
whether the prescribed maximum can be indeed reached 
by some Hamiltonian for any of the classes we have dis- 
cussed throughout this work is not difficult to answer. It 
is sufficient to take as an example the Hamiltonian induc- 
ing the Fourier transform. This transform is cyclic and its 
square leads to a permutation matrix, hence is among the 
Hamiltonians described by our method. At prescribed in- 
stances it distributes the excitation uniformly among the 
sites and hence will saturate T for any number of sites. 



V. CONCLUSIONS 

We analysed the problem of PST through passive net- 
works (channels) and the Hamiltonians which induce the 



corresponding transform. The transform is represented 
by a permutation matrix which has particularly simple 
algebraic properties. Hence, in the case of one-cycle 
permutations we were able to obtain closed expressions 
for the whole class of NN-typc Hamiltonians, while we 
proved that there exists no PST Hamiltonian in the 
framework of NN interactions. In the case of many-cycle 
permutations, the situation is more involved and for the 
time being the derivation of analytic expressions seems to 
by possible only in the case of relatively small networks. 
However, at any rate the problem can be still treated 
numerically. 

In contrast to previous work in the field, our theory is 
not limited by any a priory restrictions to the topology 
of the system and the configuration of couplings between 
different sites of the network. Hence, our approach pro- 
vides new ways for quantum wire engineering for systems 
of arbitrary topology and interactions. In this spirit, we 
were able on the one hand to recover all the known types 
of PST Hamiltonians restricted to centrosymmetric net- 
works and/or NN interaction, and on the other hand to 
derive new Hamiltonians. In principle there are infinitely 
many PST Hamiltonians which can be described within 
our theoretical framework. However, we would like to 
emphasise that the different types of Hamiltonians need 
not be always implementable by a specific experimental 
setup. Actually, which types of Hamiltonians are indeed 
implementable has to be decided on the basis of the types 
of the particular physical interactions used in the system. 
Finally, apart from studying the form of PST Hamiltoni- 
ans we also discussed the dynamics induced by them in 
terms of occupation probability distribution and overall 
bipartite entanglement. 

Additional interesting questions linked to the present 
problem such as, stability of the Hamiltonians with re- 
spect to coupling constant perturbation, detailed analysis 
of the structure of the permutation cycles and the form 
of corresponding Hamiltonians, generalisations to higher 
dimensions etc, are left for further investigation. 
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